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Lag- luminosity relation in gamma-ray burst X-ray flares: a 
direct link to the prompt emission 
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ABSTRACT 

The temporal and spectral analysis of 9 bright X-ray flares out of a sample of 113 
flares observed by Swift reveals that the flare phenomenology is strictly analogous to 
the prompt 7-ray emission: high energy flare profiles rise faster, decay faster and peak 
before the low energy emission. However, flares and prompt pulses differ in one crucial 
aspect: flares evolve with time. As time proceeds flares become wider, with larger peak 
lag, lower luminosities and softer emission. The flare spectral peak energy E p> { evolves 
to lower values following an exponential decay which tracks the decay of the flare flux. 
The two flares with best statistics show higher than expected isotropic energy E- lso 
and peak luminosity L p ^ so when compared to the E Pt i — Ei SO and E p ^ — Lj so prompt 
correlations. E Pj i is found to correlate with L lso within single flares, giving rise to a 
time resolved E p j(t) — L iso (t). Like prompt pulses, flares define a lag-luminosity rela- 



tion: L 



0.3-iokoV 



oc t . 



-0.95±0.23 



The lag-luminosity is proven to be a fundamental law 



extending ~5 decades in time and ~5 in energy. Moreover, this is direct evidence that 
GRB X-ray flares and prompt gamma-ray pulses are produced by the same mechanism. 
Finally we establish a flare- afterglow morphology connection: flares are preferentially 
detected superimposed to one-break or canonical X-ray afterglows. 

Key words: gamma-ray: bursts - radiation mechanism: non-thermal -X-rays 



1 INTRODUCTION 

The high temporal variability was one of the first 
properties to be attributed to the Gamma-ray burst 
(GRB) prompt emission in t he 7-ray energy band 
Klebsadel, Strong fc Olsonl 1 19731 ). The advent of Swift 



Gehrels et al.ll2004l ) revealed that a highly variable emission 



characterises also the early time X-ray afterglows in the form 
of erratic flares. This established the temporal variability as 
one of the key features in interpreting the GRB phenomena. 

GRB050502B and the X-ray flash 050406 
(Falcone ct al. 1 120061 : iRomano et ail [2006b ; iBurrows et al.l 
l2005bl ) provided the first examples of dramatic flaring 
activity superimposed to a smooth decay: in particular, 
GRB 050502B demonstrated that flares can be considerably 
energetic, with a 0.3-10 keV energy release comparable 
to the observed prompt fluence in the 15-150 keV band. 



Thanks to the rapid re-pointing Swift capability, it was 
later shown that flares are a common feature of the 
early X-ray aftergl ows, being present in the ~ 33% of 
X-ray light-curves l|Chincarini et al.l 120071 . hereafter C07; 
I Falcone et all 120071 . hereafter F07). On the contrary, a 
convincing optical flare, counterpart to a detected X-ray 
flare is still lacking, suggesting that the detected optical 
afterglow contemporaneous to the high-energy flares is 
dominated by a different emis sion compo nent (see e.g. 
GRB 060904B, iKlotz et al.ll2008l but see also iGreiner et all 
120091 where an optical flare was probably detected but, 
unfortunately, contemporaneous X-ray coverage is lacking). 

Based on the temporal and spectral study of a statistical 
sample of X-ray flares within GRBs, C07 and F07 showed 
that the flares share common properties and that the flare 
phenomenology can be described using averaged properties 
(see C07 and F07 and references therein): 
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• The same GRB can show multiple flares (see e.g. 
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GRB051117A which contains a minimum of 11 structures 
in the first 1 ks of observation); 

• The underlying continuum is consistent with having the 
same slope before and after the flare, suggesting that flares 
constitute a separate component in addition to the observed 
continuum; 

• Each flare determines a flux enhancement evaluated at 
the peak time AF/F between ~ 1 and ~ 1000, with a fluence 
that competes in some cases (e.g. GRB 050502B) with the 
prompt 7-ray fluence. The average flare fluence is ~ 10% 
the 15-150 keV prompt fluence; 

• Flares are sharp structures, with At/t ~ 0.1, a fast rise 
and a slower decay; 

• Each flare determines a hardening during the rise time 
and a softening during the decay time (F07 ), reminiscent of 
the prompt emission (e.g. iFord et~al]|l995h : the result is a 
hardness ratio curve that mimi cs the flare profile (see e.g. 
GRB 051117A, lGoad et alj2007t their figure 9). In this sense 
flares are spectrally harder than the underlying continuum; 

• The spectrum of a consistent frac tion of flares is bet- 
ter fitted by a Band (|Band et al.l 19931) model, simila rly to 
prompt emission pulses (see e.g. Kaneko et alj 120061. The 
flare spectral peak energy is likely to be in the soft X-ray 
range (a few keV). The spectrum evolves with time as tes- 
tified by the hardness ratio curve and by accurate spectral 
modelling. During the decay time a clear soften ing is de- 
tected (e.g. lKrimm et al.ll2007i ; iGodet et alj|2007ll ; 

• There is no correlation between the number of flares 
and the number of prompt emission pulses; 

• The vast majority of flares are concentrated in the first 1 
ks after trigger. However, late-time flares (t pca k ~ 10 5 — 10 6 
s) are present as well: flares are not confined to the steep 
decay phase, but can happen during the plateau and the 
normal decay phases. Their te mporal properties ar e consis- 
tent with those of early flares (|Curran et al.ll2008l l. even if 
their lower brightness prevents a detailed comparison with 
the entire set of early time flare properties (this is especially 
true from the spectral point of view); 

• Flares happen both in low-z and high -z environments: 
the record holder GRB 090423 at z~ 8.2 (jSalvaterra et all 
2009; lTanvir et al.| [2009l shows a prominent flare with stan- 
d ard properties when com pared to the sample of X-ray flares 
of IChincarini et all (l2010l l (C10, hereafter); 

• Flares have been detected both in hard and soft events 
such as X-Ray Flashes (e.g. XRF 050406); 

• Variability has also been detected in the X-ray after- 
glows of short GRBs (GRB with a pr ompt emission dura- 
tion Tgo < 2 s. lKouveliotou et al.J [l993l. However, given the 
lower brightness associated to these events it is still unclear 
if what is currently identified as a short GRB flare emis- 
sion, quantitatively shares the very same properties as the 
population of flares dete cted in long GRBs. GRB 050724 
l|Barthelmv et al.l [2005b) constitutes a good example of 
short GRB with late-time variability. 

From the systematic study of 113 flares in the XRT 0.3- 
10 keV energy band, as well as in 4 sub-energy bands, C10 
demonstrated that: 

• Flares are asymmetric with an average asymmetry pa- 
rameter similar to the prompt emission value; no flare is 
found rising slower than decaying; 

• The flare width evolves linearly with time w oc 0.2 t pca k. 



This is a key point which clearly distinguishes the flares from 
the prompt emission, where no e volution of the pulse-width 
has ever been found (see e.g. iRamirez-Ruiz fc Fenimorej 
E00O); 

• The width evolution is the result of the linear evolution 
of both the rise and the decay times: t r oc 0.06i pca k; id ex 
0.14t pcak . 

• The asymmetry does not evolve with time. Instead the 
rise over decay time ratio is constant with time, implying 
that both time scales are stretched of the same factor. Fur- 
thermore id ~ 2t r . Flares are self-similar in time. 

• At high energy the flares are sharper with shorter 
duration: w oc _E -0 ' 5 . Prompt pulses share the same 
property, with a sim i lar dependence on the energy band 
l|Fenimore et al.lll995l ; lNorris et al.lll996l 1; 

• The flare peak luminosity decreases with time. Ac- 
counting for the sample variance the best fit relation 
reads: L poa k oc i pc 2 a k ±0 ' 5 - The average flare luminos- 
ity declines as a power-law in time < L >oc i -1 ' 5 
l|Lazzati. Perna fc B egclman 2008); 

• The isotropic 0.3-10 keV flare energy distribution is a 
log normal peaked at ~ 10 51 erg (this can be viewed as 
the typical flare isotropic energy). The distribution further 
shows a hint of bimodality. 

• In multiple-flare GRBs, the flares follow a softening 
trend which causes later time flares to be softer and softer: 
the emitting mechanism keeps track of the previous episodes 
of emission. 

Starting from these pieces of evidence, C07, F07 and 
C10 concluded that in spite of the softness and width evo- 
lution, the X-ray flares and the prompt pulses are likely 
to share a common origin. However important questions are 
still to be addressed: do flares follow the entire set of tempo- 
ral and spectral relations found from the analysis of prompt 
emission pulses? In particular: is it possible to define a flare 
peak-lag? Do flares follow a lag-luminosity relation? Is the 
lag-time linked to other temporal properties (e.g. the flare 
duration, asymmetry, etc.)? W hat can be said about th e 
pulse start conjecture for flares (Hakkila & Nemiroff 2009)? 
Is it possible to quantify the evolution of the flare tempo- 
ral properties as a function of the energy band? Do the rise 
and decay times evolve differently with energy band? Is it 
possible to track and quantify the evolution of the spectral 
peak energy E p during the flare emission? Is there any con- 
nection between the temporal and spectral properties of the 
flares? What is the position and the track of the flares in the 
E p — Li so (spectral peak energy-isotropic luminosity) plane? 
Is there any link between the flares and the underlying X-ray 
afterglow morphology? 

This set of still open questions constitutes the major 
motivation for undertaking the present investigation. The 
primary goal of this paper is to model the X-ray flare pro- 
files and constrain their evolution with energy to obtain pa- 
rameters that uniquely qualify the shape and spectrum of 
the flares and compare those values to the well known signa- 
tures of the prompt emission pulses. Since this is often diffi- 
cult because of low statistics, overlap or non trivial estimate 
of the underlying continuum, the present work concentrates 
on very bright and isolated flares for which the underlying 
continuum does not play a major role neither from the tem- 
poral point of view nor from the spectral point of view. This 
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work is observationally driven: a critical review of theoretical 
models in the light of the present results is in preparation. 

This work is organised as follows: the sample selection 
and data reduction are described in Sect. [2] In Sect. [3] we 
perform the flare temporal analysis, while to the spectral 
properties of the sample is dedicated Sect. 3] The discussion 
follows (Sect. [5]). Conclusions are drawn in Sect. [6] 

The phenomenology of the different GRBs is presented 
in the observer frame unless otherwise stated. The 0.3-10 
keV energy band is adopted unless specified. The zero 
time is assumed to be the trigger time. The convention 
F(y,t) oc v~H- a is followed, where f3 is the spectral index, 
related to the photon index F by F = /3 + 1. All the quoted 
uncertainties are given at 68% confidence level (c.L): a 
warning is added if it is not the case. Standard cosmological 
quantities have been adopted: Ho = 70Kms -1 Mpc _ , 
Qa = 0.7, fi M = 0.3. 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

We select the brighte st, isolated flares det ected by the X- 
Ray Telescope (XRT , iBurrows et al.ll2005l ) on board Swift 
l|Gehrels et al.ll2004 ) in the time period April 2005 - Jan- 
uary 2010. The sample comprises 9 flares detected in 8 
different GRBiS GRB 050822, GRB 060418, GRB 060526, 
GRB060904B, GRB 060929, GRB070520B, GRB 070704 
and GRBO9O621A0 (see Table EU- The required isolation 
allows us to constrain the flares temporal and spectral prop- 
erties with high accuracy; their brightness assures the pos- 
sibility to analyse in detail the evolution of their temporal 
properties in different energy bands inside the 0.3-10 keV of 
the XRT; at the same time, this requirement guarantees that 
the partial ignorance of the underlying continuum proper- 
ties has a negligible impact on our conclusions, both from 
the spectral and from the temporal point of view. 

XRT data have been processed with the HEASOFT pack- 
age v. 6.6.1 and corresponding calibration files: standard 
filtering and screening criteria have been applied. Swift- 
XRT is designed to acquire data using different observing 
modes to minimise the presence of pile-up. When the source 
is brighter than a few counts -1 , the data are acquired in 
Windowed Timing (WT) mode. In this case we applied 
th e standard pile - up corr ections following the prescriptions 
ofH omano et all i|2006ah when necessary. For lower count- 
rates the spacecraft automatically switches to the Photon 
Counting (PC) mode to follow the fading of the source. The 
events are extracted from circular regions centred at the af- 
terglow position with progressively smaller radii to assure 
the best signal to noise ratio (SN). When the PC data suf- 
fered from pile-up, we extracted the source events in an an- 
nulus whose inner radius is derived comparin g the observed 
to the nominal point spread function (PSF. iMoretti et alj 

1 GRB 050502B is not included since the bright flare detected by 
the XRT is possibly due to the superposition of more than one 
structure iBurrows et al.ll2005bl ; C10) 

2 Note that for GRB0 90621A, the Burst Alert Tel escope (BAT 
Barthe lmv et al. 2005a) triggered on the precursor l lCurran et al.l 
20091) . 
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Figure 1. Best fit profile of the flare detected in GRB060904B in 
different XRT energy bands. The thick line, dotted line, dashed 
line and dot-dashed line refer to the flare profile observed in the 
0.3-1 keV, 1-2 keV, 2-3 keV and 3-10 keV energy band, respec- 
tively. The profiles have been re-normalised for the sake of clarity. 
The blue arrows point to the flare peak times in the different en- 
ergy bands. Softer energy bands profiles peak later. Inset: The re- 
normalised best fit profiles have been aligned at tpeak = 170.71 
s, peak time of the flare in the 0.3-10 keV energy band as de- 
rived from the best fit parameters of Table IA1I Softest profiles 
are wider. 

120051 : IVaughan et al.ll200rj ). The background is estimated 
from a source-free portion of the sky and then subtracted. 
The background subtracted, PSF and vignetting corrected 
light-curves were then re-binned so as to assure a mini- 
mum signal-to-noise (SN) equals to 4 for WT and PC mode 
data. Data coming from different orbits of observation were 
merged to build a unique data point, when necessary. This 
procedure was applied to extract GRB light-curves in the 
nominal XRT energy range (0.3-10 keV) as well as count- 
rate light-curves of the same event in 4 different sub-energy 
bands: 0.3-1 keV; 1-2 keV; 2-3 keV and 3-10 keV0. 



3 FLARE TEMPORAL ANALYSIS 

Each flare profile in each energy band is modelled using the 
iNorris et all (120051 ) (hereafter N05) profile. This choice al- 
lows a direct comparison between the properties of the X-ray 
flares and of the wide, long-lag prompt pulses analysed by 
N05 (see Sec.©. The N05 profile reads: 

I(t) = A\exp[-— ^— - for t > t s (1) 

(I — I s ) T2 

where A = exp(2/x) and fi = (ti/t2) 1//2 . t 3 is the pulse onset 
time. The parameters are defined following N05. In partic- 
ular, the intensity peaks at: 

3 Note that, differently from C10 the described data reduction 
assures a minimum SN equal to 4 for each data point of the GRB 
light-curves in the different energy bands. This allows a higher 
sensitivity which translates into a more accurate evaluation of 
the best fit parameters. 
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Figure 2. Panels (a) through (e): GRB 060904B flare best fit profile in the total 0.3-10 keV XRT energy band and in the 4 channels. 
Insets: BAT signal contemporaneous to the flare emission detected in the 15-150 keV, 15-25 keV, 25-50 keV, 50-100 keV and 100-150 
keV energy bands. The typical uncertainty affecting the BAT data is also shown. Panels (f) through (m): observed evolution of the 
flare decay time (f), rise time (g), width (h), asymmetry (i), start-time (I) and peak time (m) with observed energy band. The best fit 
power-law relation describing the evolution of each parameter is drawn with a red solid line and explicitly written in each panel. The fit 
rpnnrtpd in nanrl fm ) concerns the XRT data onlv 
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£peak = t s + (TlT2) 1/2 (2) 

The pulse width is measured between the two 1/e intensity 
points and consequently reads: 

™ = T 2 (1 + 4^) 1/2 =t r +t d (3) 

The asymmetry is denned as: 

k = (i+^r 1/2 = r^- ( 4) 

Id + tr 

where td and t r are the 1/e decaying and rising times, re- 
spectively: 

td, r = ir 2 [(l + 4 M ) 1 / 2 ±l] (5) 

We account for the entire covariance matrix during the er- 
ror propagation procedure. The continuum underlying the 
flare emission is estimated from data acquired before and af- 
ter the flare and is modelled using a power-law or a broken 
power-law: the parameters related to the continuum have 
been left free to vary in a first fit and then frozen to their 
best fit values in a second fit. The best fit parameters and re- 
lated quantities are reported in Table IAT1 Figure [1] shows the 
best fit profiles in the different XRT band-passes obtained 
for the flare detected in GRB 060904B, taken as an example. 
A general trend can be seen for the flare to be wider and 
peak later at lower energies. A complete summary of the re- 
sults obtained from the fitting procedure in different energy 
bands is portrayed in Fig. [2] the detected evolution of the 
different parameters with energy (panels / to m) will be ex- 
tensively treated in SubSec. 13.21 and 13.31 The flare detected 
in GRB 060904B is here shown as an example. 

3.1 Best fit parameters correlations 

Figures [3] and [4] show the correlations or the lack thereof, 
for the quantities derived from the best fit parameters of 
Table IA1I superimposed to the results obtained by N05 for 
a sample of long-lag wide prompt GRB pulses and by CIO 
for a sample of 113 early X-ray flares. In these plots we 
see that bright flares obey the set of correlations found for 
the complete X-ray flare catalogue: the rise time is linearly 
correlated with the decay time, and both linearly evolve with 
time. This gives rise to a flare width that linearly grows with 
time (see Fig. [3j . 

Figure [4] depicts instead the relation between the asym- 
metry k and the temporal parameters describing the flare 
profiles and demonstrates that it is not possible to distin- 
guish the bright from the main sample on the basis of their 
asymmetry. Furthermore, the asymmetry seems to be inde- 
pendent of the other parameters (Fig. |4]). 

3.2 Best At parameters evolution with energy 

The well known trend of narrower pulses at higher energy 
is apparent in Fig. [1] for the flare detected in GRB 060904B. 
This trend is quantified in Fig. [5] for the entire sample: 
the narrowing with energy follows a power-law behavioui0: 



w oc E °' 32 ± ' 01 The same is true for the rise and decay 
times: t r cx E - - 37±OM ; td oc S" - 39 * - 01 . The width evo- 
lution is the combined result of the evolution of these two 
time scales and its evolution is consequently associated to 
a power-law index which is between 0.29 and 0.37. Notably 
the rise time shows a steeper dependence on the energy band 
than the decay time. As a result, the asymmetry is slightly 
dependent on energy: flares are on average more asymmetric 
at higher energies. Flare profiles peak later at lower energies 
(Fig. [2] panel m); on the contrary, it seems that high energy 
profiles start their rise later (Fig. [2] panel I): while the flare 
in GRB 060904B is here portrayed as an example, the same 
is true for the entire sample as testified by the results of 
Tabled] 

Figure [S] also demonstrates that the energy dependence 
of the various parameters varies from one flare to another: 
the average relations describe the general behaviour of the 
parameters with energy band, but the evolution in a par- 
ticular flare can be markedly different. This statement is 
quantified in Table [T] where the best fit power-law index 
describing the behaviour of the temporal parameters with 
energy band is listed. The best fit relations are portrayed in 
Fig. |1 panels (f) through (m) for GRB060904B. How this 
different evolution rate relates to the spectral properties of 
the flares is investigated in Sec. [4] 

3.3 Pulse peak lag 

A commonly used parameter related to the spectral evolu- 
tion of the pulses is the time-lag between two energy chan- 
nels. In this work we always refer to the flare peak lag which 
is defined as the difference between the peak of the lowest 
and highest energy profiles measured by the XRT as derived 
from the best fit parameters of Table IA1I Taking the flare 
detected in GRB 060904B as an example, from Fig.[T]a clear 
tendency is apparent for the peak in the harder channel to 
lead that in the softer channel: this translates into a measur- 
able positive pulse peak lag. The same is true for the entire 
sample of flares. The evolution of the flare peak times as a 
function of energy in the XRT bandpass is reasonably well 
represented by a power-law, as testified by Table [T] and by 
Fig. [5] Again, different flares have different rates of evolution 
of their peak times in fixed energy band-passes. 

We scanned the BAT data looking for detections of the 
bright flares in the 15-150 keV energy range: a peak finding 
algorithm is run. The software is written to automatically 
scan the observations using all the possible re-binning time 
scales: the central time of the bin which maximises the de- 
tection is used as flare peak time. The flare peak times in 
the 7-ray band have been added to Fig. [2] panel (m) when 
a 5cr detection was found. As appears from this figure, the 
flare peaks at high energy are not consistent with the ex- 
trapolation of the law deduced from X-ray data alone. The 
same is true for the flares in GRB 060418, GRB 060526 and 
GRB 070704 which have a BAT detection contemporaneous 
to the flare emission. With reference to GRB060904B, it is 
notable that the highest rate of evolution of i pea k is found for 



4 For the 0.3-1 keV, 1-2 keV, 2-3 keV energy bands, we define an 

/J* 1 Ef(E)RM(E)d{E) 

where RM stands 



effective energy E e g 



Ie2 f(E)RM(E)d(E) 



for the XRT response matrix, f(E) ~ E~ 1 , as found by C10 
for the average flare spectral energy distribution. For the hardest 
band the best fit spectrum of each flare is used. 
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Figure 3. Scatter plot of the temporal parameters describing the flare profile as derived from the best fit values reported in Table fATl 
(a): decay vs. rise time for our sample of bright flares (black dots), for the sample of flares analysed by CIO (grey dots) compared to the 
values found by N05 for a sample of long-lag, wide prompt pulses detected by BATSE. Blue dashed line: locus of the points for which 
*r = td- Red solid line: best fit law for the three samples: t<j = (1.1 ± 0.1) + (1.9 ± 0.1)t T . (b): decay vs. peak time for the three samples 
described for panel (a). Red solid line: linear best fit relation found for the flare samples: = (3.0 ± 0.4) + (0.14 ± 0.02)t pca ]j. (c): rise 
vs. peak time. The best fit relation found for flares is marked with a red solid line: t T = (4.0 ± 0.3) + (0.06 ± 0.01)t pca ] I . (d): width vs. 
peak time. Red solid line: flare best fit relation: w = (4.8 ± 0.6) + (0.20 ± 0.02)t pea k. 



Table 1. Evolution of the flare profile with energy band. The rise time, decay time, peak time, width and asymmetry have been fitted 
using a power-law model in energy: y oc E a . From left to right: power-law index associated to the evolution of the following parameters: 
rise time at r , decay time ce td , peak time a tp , width ct w and asymmetry a^. 



GRB 




atd 






a w 


Oik 


050822 


-0.85 ±0.13 


-0 


56 ± 0.07 


-0.11 ±0.03 


-0.83 ±0.10 


0.12 ±0.04 


060418 


-0.60 ±0.04 


-0 


69 ± 0.02 


-0.021 ± 0.005 


-0.68 ±0.02 


-0.034 ± 0.009 


060526f^ 


-0.46 ±0.02 


-0 


43 ±0.03 


-0.029 ±0.001 


-0.44 ±0.02 


0.040 ± 0.005 


060526s^ 


-0.73 ±0.08 


-0 


40 ±0.02 


-0.042 ± 0.006 


-0.43 ±0.01 


0.24 ±0.01 


060904~ 


-0.41 ±0.02 


-0 


40 ±0.02 


-0.090 ±0.014 


-0.41 ±0.06 


0.0060 ± 0.039 


060929 


-0.21 ±0.02 


-0 


17 ±0.02 


-0.050 ± 0.006 


-0.20 ±0.04 


0.090 ±0.018 


070520B 


-0.67 ±0.06 


-0 


32 ±0.03 


-0.063 ± 0.025 


-0.41 ±0.07 


0.48 ±0.05 


070704 


-0.21 ±0.04 


-0 


21 ±0.03 


-0.033 ±0.011 


-0.18 ±0.04 


0.039 ±0.021 


090621A 


-0.19 ±0.04 


-0 


21 ±0.02 


-0.029 ±0.021 


-0.21 ±0.18 


-0.028 ±0.057 



a First flare detected by the XRT; 
6 Second flare detected by the XRT. 
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Figure 4. Asymmetry vs. rise time (panel (a)), decay time (panel (b)), width (panel (cj) and peak time (panel (d)) for our sample of 
bright X-ray flares (black dots), the sample of flares from C10 (grey dots) and for a sample of long-lag, wide prompt pulses from N05. 
In each plot, the arrow points to the direction of growing asymmetry. 



band-passes which are crossed by the spectral peak energy 
E p during the flare emission (see SubSec. l4.3p . 

The presence of correlations between the peak-lag and 
the flare shape parameters is investigated in Fig. [7] The 
lag is found to be positively correlated with both the rise 
and the decay times (panel a), giving rise to the width- lag 
correlation portrayed in panel b: the wider the flare, the 
higher the lag value. The same is true for the flare peak 
time: later flares are associated to larger time lags. 



4 FLARE SPECTRAL ANALYSIS 

The flare emission has been proven to undergo a strong spec- 
tral evolutio n with the hardne ss ratio tracking the flare pro- 
file (see e.g. lGoad et ai]|2007r ). This results in flares profiles 
on average harder than the underlying continuum. In this 
section we quantify the spectral properties of each flare of 
our sample using the hardness ratio analysis (Subsec. 14. 1[) . 
A link between the temporal and spectral properties is es- 
tablished in Subsec. 14.21 while a proper spectral modelling 
is performed for the two flares with the best statistics, mea- 
sured redshift and BAT detection in Subsec. 14.31 



4.1 Hardness ratio analysis 

We define the total hardness ratio as the ratio between the 
re-constructed counts in the 0.3-2 keV and 2-10 keV energy 
bands, calculated in the time interval t\ through tf, where 
ti = tpcak —ti and if = t pca k + td ■ £ pca k, t T and t d are derived 
from the best fit parameters of the 0.3-10 keV profile: 



HR t . 



Counts(2 - lOkeV) 



Counts(0.3 - 2keV) 



(6) 



During the rise and the decay of each flare, the hardness 
ratio is defined as follows: 



HRrisc — 



Counts(2 - 10 keV) 



Counts(0.3 - 2keV) 



HR 



decay 



Counts(2 - lOkeV) 



Counts(0.3 - 2keV) 



(7) 



(8) 



Figure [8] clearly shows that flares are harder during the 
rise time and softer during the decay time. The inset of 
this figure illustrates the presence of a trend: the harder the 
overall profile, the lower is the spectral ratio between the 
rise and the decay portions of a flare. A lower degree of 
spectral evolution is detected for spectrally harder flares, in 
agreement with the findings of Subsect. 14.21 Note however 
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Figure 7. (a): Bright flares rise (red triangles) and decay times (black dots) vs. peak lag. The red dot-dashed line and black dashed 
line indicate the best fit relations: t r = l(r a49±0 ' 28 (tlag) 1 - 3 * ' 21 ; t d = lO" - 31 * ' 32 (ti ag ) 15±0 - 24 . (b): Width vs. peak lag. The dashed 
line indicates the best fit relation: w = 10 ' 15±0 ' 27 (tiag) 1 ' 2 ^ ' 20 . (c): peak time vs. peak lag. The dashed line indicates the best fit 
relation: tpcak = 10 ' 69 (ii a g) 1,2±C1 ' 21 . (d): scatter plot of the asymmetry vs. peak lag. No relation is apparent. The arrow points to 
the direction of growing asymmetry. 



that the difference between the HR r i sc and HRdccay is higher 
for harder HR tot . 

The connection between the temporal and spectral 
properties is made in Fig. [9] Panel (c) clearly shows that 
hard flares are less asymmetric. 



4.2 Rate of evolution of a flare profile with energy 
vs. HR 

The inset of Fig. [S] suggests the presence of a correlation 
between the flare observed spectral evolution and the over- 
all spectral hardness: harder flares shows a limited level of 
detected spectral evolution. In this subsection we prove that 
a similar conclusion is reached through the analysis of the 
rate of evolution of the temporal parameters describing the 
flare profiles. 

The evolution with energy of each flare profile has been 
analysed in Subsect. l3.2| parametrising the evolution of each 
flare-shape parameter (rise time, decay time, peak time, 
width and asymmetry) using a power-law model. The best 
fit power-law indices are reported in Table[T]and used as evo- 
lution indicators in Fig llOl This figure clearly shows that the 
rate of variation of a flare temporal profile is a function of 
the detected hardness of the emission: the harder the emis- 
sion, the lower is the evolution. The parameters describing 



the profile of softer flares undergo a substantial evolution 
from one energy band to the other (corresponding to higher 
\a\ values), while the evolution of hard profiles is limited 
(lower | a | values). Figure 1101 panel (b) confirms that the 
evolution of the rise time from one energy band to the other 
is more sensitive than the decay time to the average spec- 
tral hardness of the flare. Furthermore, it can be seen that 
t T and td undergo a similar evolution with bandpass only in 
hard flares: the softer the flare, the higher is the difference 
between the i r and t<j rate of evolution with energy, with 
the rise time always evolving faster in energy. The presence 
of the trends depicted in the first and second panels from 
the top of Fig. 1 101 automatically translates into the results of 
the third one: since w — t r + id, it is not surprising that the 
width is subject to a higher level of evolution in softer flares. 
No clear correlation is found between the rate of evolution 
of the asymmetry and the spectral hardness: however, it is 
worth noting that the asymmetry evolves in the four softer 
flares, while a limited level of evolution or no evolution is 
detected in the five hardest flares (shaded area of Fig. 1101 
bottom panel). This is a consequence of the results above. 

At this point, a lag-hardness correlation is expected: 
Fig. [11] demonstrates that softer flares are characterised by 
larger time lags. 
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Figure 9. Best fit observed width (a), peak time (b) and asymmetry (c) as a function of the total hardness ratio. The red arrow in 
panel (c) points to the direction of higher asymmetry. 



Table 2. Average spectral properties of the flares detected in 
GRB 060418 and GRB060904B. The two spectra have been ex- 
tracted in the time interval defined by (t pca k — ir)-(< P cak + *d) 
which corresponds to the shaded area in the insets of Fig. 1121 
The fitting model is a photo-electrically absorbed Band func- 
tion with E p as free parameter (tbabs*ztbabs*ngrbep within 
xspec). The intrinsic hydrogen column density TVjj z is frozen to 
the value found from the joint fit of the time resolved spectra (see 
Table lA2l and lA3l . The Galactic absorption in the direction of the 
bursts has been accounted for following Kalberla et al., 2005. The 
isotropic equivalent energy i?i so , the isotropic equivalent luminos- 
ity L; ao and the isotropic equivalent peak luminosity T Pi i so are 
computed in the 1 — 10 4 keV rest frame energy range. Z/ p> i so is 
calculated re-scaling the time integrated spectrum with the flare 
peak flux. Errors are provided at 90% c.l. 





GRB 060418 


GRB 060904B 


z 

Pb 

E p (keV) 

Afa.z (10 22 cm" 2 ) 

E iso (erg) 

L iso (erg s _1 ) 

£ p ,iso(erg s _1 ) 


1.489 
-0.9 
-2.30 ±0.09 

n n + 0.2 

°-°-0.3 

0.59 

(5.4 ±0.5) x 10 51 
(4.1 ±0.4) x 10 50 
(5.1 ±0.5) x 10 50 


0.703 
-1.0 
-2.47 ±0.10 

o 7+O.2 
°- ' -0.2 

0.52 

(2.4 ±0.2) X 10 51 
(5.2 ±0.5) x 10 49 
(9.1 ±0.9) x 10 49 


X 2 /dof 
P-value 


209.26/216 
0.61 


312.86/304 
0.35 



4.3 Spectral modelling 

The flare emission is typical of the X-ray regime. The lim- 
ited energy window of the XRT (0.3-10 keV), as well as the 
complete ignorance of the distance of the emitting source 
coupled to a poor knowledge of the intrinsic neutral hy- 
drogen column density absorbing the softest photons, can 
lead to the misidentification of intrinsically curved spectra 
in terms of absorbed power-laws. We therefore select a sub- 
sample of flares with measured redshift and with contem- 
poraneous detection by the BAT in the 15-150 keV: these 
requirements restrict the analysis to the two flares observed 
in GRB060904B and GRB 0604183. 

The very good statistics characterising both signals 



5 In spite of the measured redshift and BAT detection, 
GRB 060526 is excluded because of the partial overlap of the two 
structures. 



gives us the possibility to perform a time resolved spec- 
tral analysis. To follow the spectral evolution of the source, 
we time sliced the BAT and XRT data of GRB 060418 and 
GRB 060904B into 4 and 8 bins, covering the 129-320 s and 
122-170 s time intervals, respectively. The time intervals 
were defined so as to contain a minimum of ~1000 pho- 
tons in the XRT range and to possibly have enough signal 
in the BAT energy range to constrain the spectral parame- 
ters. The spectra were first fitted using an absorbed simple 
power-law (SPL) model within XSPEC: the model is absorbed 
by both a Galactic and an intrinsic absorption components. 
The Galactic compon ent was frozen to the value specified by 
iKalberla et al.l (|2005t ) in the direction of the burst (which is 
taken to be the UVOT-XRT refined position delivered by 
the Suii/i-XRT team); the intrinsic neutral hydrogen col- 
umn density A?h,z was left free to vary. When possible we 
took advantage of the simultaneous BAT and XRT coverage, 
performing a joint BAT-XRT spectral fit. The normalisation 
for each instrument was always tied to the same value. The 
results are reported in Table [M] and Table El The SPL fit 
results are statistically unsatisfactory and are not able to ac- 
count for the 15-150 keV data. Moreover, the best fit A?h,z 
are inconsistent with the values obtained for both GRBs 
from late time spectra, where no spectral evolution is de- 
tected and no spectral break is expected to lie in the XRT 
band pass: for GRB 060904B and GRB 060418, from spec- 
tra extracted in the time intervals 1-500 ks and 200-1000 
s, respectively, we obtained iVa> ~ (0.4 ± 0.1) x 10 22 cm -2 
(90% errors are provided). The SPL fits of Table Upland [Ml 
instead indicate values > 10 22 cm -2 . These facts, together 
with the clear evolution of the best fit spectral photon index 
to softer values with time, suggest that the peak of the vF v 
spectrum is moving through the BAT+XRT bandpass. We 
tested this possibility performin g time resolved sp ectral fits 
to the empirical Band function l|Band et al.lll993l ) using E p 
as free parameter (ngrbep model). The high and low power- 
law indices (qb and /3b) which characterise the Band model 
were us ed as free paramete rs: we froze q:b to the typical -1 
(see e.g. lKaneko et aill2006l ) value when poorly constrained. 
The uncertainties affecting each parameter were computed 
freezing aa to the best fit value. The spectra of each flare 
were fitted simultaneously under the assumption of a con- 
stant A^h.z during the flare emission. The results of the fit- 
ting procedure are listed in Table|M]andrS2]for GRB 060418 
and GRB 060904B, respectively. The observed evolution of 
the peak energy with time is represented in Fig |12l for both 
flares, Ep is found to evolve to lower values, with a decay 
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Figure 5. Evolution of the width (a), rise time (b), decay time 
(c) and asymmetry (d) with observed energy band obtained for 
our sample of bright X-ray flares. Results coming from the fit of a 
particular flare in different energy bands are indicated by dots and 
linked by a dotted line. Red stars: mean values derived from the 
entire sample. Red solid line: best fit power-law model: (a) w = 
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Figure 6. Best fit power-law index describing the evolution of 
the decay time with observed energy as a function of the best 
fit power-law index found for the evolution of the rise time with 
energy. Red dashed line: locus of the points for which a t( j = atr- 
Darker points are associated to spectrally harder flares. 
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Figure 8. Hardness ratio of each flare during the rise time 
vs. hardness ratio during the decay time: the flare spectrum is 
harder during the rise time. Red solid line: best fit linear rela- 
tion: HR risc = (0.07 ±0.01) + (1.6 ±0.1)HR docay . Inset: rise time 
over decay time hardness ratio vs. overall spectral hardness of the 
flares: a trend is apparent for harder flares to show more similar 
hardness ratios during the rise and decay portions of X-ray the 
light-curves. 



that during the flare decay time can be represented by ei- 
ther a power-law with index a ~ —7 or an exponential with 
e-folding time r = 29 ± 5 s (GRB 060904B) and r = 20 ± 2 
s (GRB 060418). The uncertainty of the inter-calibration of 
the BAT and XRT has been investigated as possible source 
of the detected spectral evolution: for each time slice, we 
multiplied the fit model by a constant factor which is frozen 
to 1 for the BAT data. For XRT, this factor is left free to 
vary between 0.9 and 1.1, conservatively allowing the XRT 
calibration to agree within 10% with the BAT calibration. 
The best fit parameters found in this way are completely 
consistent with those listed in Table lA"2l and lA"3l The inter- 
calibration is therefore unlikely to be the main source of the 
observed evolution. 

For each time interval we computed the isotropic equiv- 
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Figure 10. Best fit power-law index describing the evolution of 
the flare temporal profile parameters with observed energy (see 
Table U) as a function of the spectral hardness of the flare. The 
evolution of each flare parameter with energy band has been fit 
using the model: y oc E a . (a): rise time power-law indices, (b): 
black dots: decay time power-law indices; grey dots: rise time 
power-law indices shown for direct comparison, (c): Width power- 
law indices, (d): asymmetry power-law indices. The shaded area 
represents a region of low evolution. In each panel an arrow points 
to the direction of stronger spectral evolution of the parameter. 
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Figure 11. Observed pulse peak lag as a function of the spectral 
hardness of the flares. A power-law with decaying index -0.5 is 
also shown for comparison (red dashed line). 



alent luminosity Li so in the rest frame 1-10 4 keV energy 
band: these contributions have been de-absorbed for both 
the Galactic and intrinsic hydrogen column density (see Ta- 
ble IA2I I A3f) . The spectral peak energy is found to evolve 
in the Li so — E P: i plane following a power-law behaviour 
with best fit power-law index ~ 1 (Fig. 1 1 3p . In particular, 
for GRB 060418 (GRB060904B) we have E Pri oc Lf s ° ±0 - 2 
and (£ p ,i oc I/^ ' 2 )). For GRB060904B we were able to 
extract a spectrum at the onset of the flare emission: the 
properties of this spectrum (green cross in Fig. 1131 upper 
panel) are clear outliers with respect to the general trend 
defined by the other spectra. The same can not be done for 
the GRB 060418 flare: in this case, a spectrum extracted at 
the onset is completely dominated by the contribution of 
the underlying continuum and nothing reliable can be said 
about the flare spectral properties in this time interval. 

The average spectral properties of the two flares have 
been computed extracting a spectrum from i pca k — t z 
to tp Ca k + id (shaded area of the insets of Fig. I12p : in 
both cases the BAT detection rules out a simple power 
law behaviour (x 2 /dof=289.34/216 P-val=6 x 10~ 4 and 
X 2 /dof=466. 45/304 P-val=6 x 10~ 9 , for GRB 060418 and 
GRB060904B, respectively). The two spectra are instead 
well represented by a Band model with best fit observed 
peak energies E p = 3.32.Q' 2 . keV (intrinsic value £? p ,i = 8.2 
keV) for GRB 060418 and E p = 3.7lo. 2 keV (E p ,[ = 6.3 
keV) for GRB060904B (90% uncertainties are provided). 
The best fit results are reported in Table [2] note that 
the isotropic equivalent peak luminosity L P] i S o is calculated 
starting from the time averaged spectrum re-scaled with the 
flare peak flux. In spite of not representing the real flare 
peak luminosity, this definition allows us to perform a di- 
rect comparison to the results obtained in the literature 
for the gamma-ray pr ompt pulses (see e.g. iNava et~aL I l2008l . 
lYonetoku et alj|2004l ). 



5 DISCUSSION 

Through the analysis of 9 bright X-ray flares we proved the 
existence of a set of correlations which link both the tempo- 
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Figure 14. E p \ — E- lso and E Pi i — L p ^ so planes as obtained for the two X-ray flares detected in GRB 060904B and GRB 060418 (blue and 
red dot, respectively) superimposed to the data obtained from the sample of 83 GRBs with measured redshift and spectral parameters 
analysed in Nava et al, 2008. Black solid lines: best fit to the enlarged sample of 110 GRBs published in Ghirlanda et al., 2009, with best 
fit slopes 8 = 0.476 ± 0.025 and <5 = 0.395 ± 0.024 for the left and right panel, respectively. The dashed lines represent the 3<x scatter. 
In both panels a red cross marks the prompt spectral properties of GRB 060418. The approximate position of short GRBs and XRFs in 
the E p i — Ei so plane is marked with shaded areas. 



Table 3. Summary of the relations which link the flare parameters, f (4-) stands for the presence of a direct (inverse) proportionality. 
means that no correlation has been found T Pi i so is in units of ergs -1 ; temporal parameters are indicated in s. Errors are reported at 90% 
c.l. for correlations involving the peak luminosity. 
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Table 4. First and second columns: summary of the best fit re- 
lations describing the parameters evolution with energy band. 
Third and fourth columns: rate of evolution |aj of the parame- 
ters with energy as a function of the spectral hardness: a 4- (t) 
indicates inverse (direct) correlation. means that no correlation 
has been found. 
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ral and spectral flare parameters. The choice of this sample 
introduces an obvious caveat in our analysis, namely, that 
we make the explicit assumption that the results we derive 
from the bright flares hold for all the flares. However, apart 
from the brightness, there is no evidence that our sample of 
bright flares is not representative of the entire population 
at the 80% c.l. as obtained by a Kolmogorov-Smirnov (KS) 



test. We therefore conclude that our results are likely to be 
extended to the entire population of X-ray flares. 



5.1 Temporal parameters evolution with energy 

Flares have different temporal profiles in different energy 
bandpasses: in strict analogy to prompt emission pulses, 
flares are broader and peak later at lower energy. C10 found 
w oc E~ ' 5 . Our sample of bright flares shows w oc £' _0 ' 3 . 
For the prompt puls es the relati on reads: w oc E~ 4 
(IFenimore et al. 1 1 19951. but see also [Borgonovo et al. 1 120071 . 



N05, iNorris et al.l 119961 and references therein). The width 



evolution results from the joint evolution of two time scales: 
the rise and the decay times. The evolution with energy of 
tr and td is well described by a power-law, with power-law 
indices which are found to be very similar (see Table [T]): 
however, the very good statistics of our sample allows us to 
conclude that the rise time is more sensitive to the energy 
bandpass than the decay time. This is evident from Fig. HJ] 
Since the evolution of the temporal parameters with energy 
band is strictly linked to the spectral evolution during the 
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Figure 12. Temporal evolution of the observed peak energy of 
the vF v spectrum of the flare detected in GRB 060904B (panel 
(a)) and in GRB 060418 (panel (b)). Blue dashed line: best 
fit power-law decay for the evolution of E p during the decay 
time. The best fit power-law decay index is ay = —7.2 ± 0.7 
(GRB 060418) and a d = -6.5 ± 0.7 (GRB 060904B). Red dot- 
dashed line: best fit exponential decay of E p during the flare 
decay time. The zero time of the power-law is assumed to be 
the BAT trigger time. The best fit e-folding time is found to be 
r = 20 ± 2 s and r = 29 ± 5 s for GRB 060418 and GRB 060904B, 
respectively. Insets: Temporal evolution of E p compared to the 
evolution of the 0.3-10 keV flare profile (red solid line) . The flare 
flux is represented in arbitrary units. The shaded areas mark the 
interval of time of extraction of the average flare spectrum, from 
ipcak — ti to tp ea k + td- ln each panel, the orange dotted lines 
mark the flare peak time as derived from the best fit parameters. 



flare emission, this would point to some spectral differences 
between the rise and decay portions of the light-curve. In 
particular, the different relative position of the evolving E p 
with respect to the fixed energy bands during the t T and is 
likely to play a major role: this can be in turn a footprint of 
two different physical mechanisms dominating the rise and 
the decay phases. The first consequence of this finding is that 
each energy channel is only nearly an exact time-stretched 
version of the others: flares can be treated as self-similar in 
energy only at the 0th level of approximation. The second 
direct consequence is the presence of a slight trend pointing 
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Figure 13. Time resolved E p ^ — L- lso correlation for the flares 
detected in GRB060904B (upper panel) and GRB 060418 (lower 
panel) . Red solid lines: best fit power-law relations: E p ^ oc 



oc Lt: 



to more asy mmetric flar e s at h igher energy (Table [4] Fig. 
panel cO. lNorris etTafl (|l996h defining the asymmetry as 
i r /id found no evolution of this ratio with energy for GRB 
prompt pulses: however, the separate evolution t r (E) and 
td(E) was not calculated. 

While the soft energy flare profiles are clearly broader 
and peak later, it is unclear if the flares start simulta- 
neously at all ene rgies. The Pulse Start Conjecture was 
first proposed by 



Nemirofj (|2000l ) and then tested by 



lHakkila fc Nem irofi (200^): these authors studied a sample 
of 199 prompt BATSE pulses in 75 different GRBs using an 
automatic pulse fitting methodology. Each pulse is modelled 
by a N05 profile (the same profile used for flares). We caution 
that a strong coupling is built in the N05 profile between t s 
and TV. the two parameters are clearly related since both fit 
the flares before the peak. This translates in a non-negligible 
probability of obtaining a good fit with perfect \ 2 but with 
non-physical t s and n: any slight but systematic variation 
t s (E) would be destroyed. The other 2 major problems are 
connected to the structure overlap a n d low statistics that 
affect part of the lHakkila fc Nemirofj (|2009l ) sample as dis- 
cussed by the authors. In spite of these limitations the Pulse 
Start Conjecture was concluded to hold for prompt gamma- 
ray pulses with an uncertainty of ~ 0.4 s. For the X-ray 
flares the situation is different: the selection of 7 isolated 
flares with very good statistics led to the conclusion that a 
slight trend for flare profiles at high energy to start later is 
present. An example is shown in Fig. [2] panel I. Table lATI 
reports the start time values for the entire sample of bright 
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flares: the high energy profile t s is always larger than the 
low energy t s : in 5 cases out of 7 the difference between the 
two values is significant at more than 3 a c.l. 

The flare emission preferentially builds up at lower en- 
ergies at the beginning of the GRB flares. 

Figure [S] demonstrates that the rate of evolution of the 
temporal parameters with energy in a single flare can be 
markedly different from the average sample behaviour (see 
the values reported in Table[T]compared to those of Tabled}. 
Notably, this is directly linked to the observed flare spectral 
properties as proven by Fig. 1101 the harder the flare, the 
lower is the detected rate of evolution of the flare profile 
in different energy bands. Since hard flares also display a 
limited spectral evolution in the 0.3-10 keV energy range 
(insets of Fig. [SJ, (likely because their E p lies outside the 
XRT energy band during the majority of the emission time), 
the result above implies a strict link between the spectral 
properties and the evolution of the temporal properties. 

The same is true when the flare peak lag is considered: 
in the case of the prompt emission, the cross-correlation lag 
between energy bands was shown to anti-correlate with the 
BATS E spectral hardness ratio bv lNorris. Marani fc Bonnell 
( 2000) ; when pulses properties are considered instead of time 
integrated burst properties, the result is that the pulse spec- 
tral hardness anti-correlates with pu lse lag and duration , 
and correlates with pulse intensity (jHakkila et al.l [2008). 
Figure [11] demonstrates that X-ray flares share the same 
property: a trend is evident for flares with larger HR to t to 
have shorter peak lags. The lag-hardness correlation of this 
figure, together with Fig.|8]again implies a robust connection 
between the flares temporal properties and their spectral 
evolution: a stronger spectral evolution directly translates 
into a higher lag. This means that the lag characterising 
the flares profiles is the direct result of t heir spectral evolu- 
tion. T he same conclusion was reached bv lKocevski fc Liand 
(|2003h for the prompt emission: these authors concluded 
that the fundamental origin of the observed lag is the evo- 
lution of the GRB prompt spectra to lower energies. Our 
results of Fig. \W\ extend this finding to other flare temporal 
properties: the rate of evolution with energy of the rise time, 
decay time and width is inversely correlated to the spectral 
hardness and directly linked to the spectral evolution. 

Finally we mention that only soft flares have been 
found to be highly asymmetric (Fig. [2] panel c). Similarly 
iNorris et al.l l| 19961 ) found that soft prompt pulses are on 
average more asymmetric. 



5.2 The flare lag-luminosity relation 

The flare peak lag is directly related to the flare duration 
(Table|3] Fig.O and since w has been found to vary inversely 
with the flare peak luminosity (C10, but see also the relation 
reported in Table [3} , an inverse correlation between peak 
lag and luminosity is expected. Since only three flares in our 
sample have measured redshift, we selected the flares with 
the best statistics, measured redshift and well defined peak 
times in both the hardest and softest XRT energy bands 
from the sample of C10. Figure [T51 confirms the expectation 
above, showing that the X-ray flares define a rather tight 
lag-luminosity relation: 



Table 5. Prompt 7-ray lag-luminosity relations. 



Sample Rcf. 



L P,i = o _ iq(51. 5±0.4) ( 'l»8. 1 
(ergs- 1 ) ^ s 

L p,iso _ in(51-54±0.05) I 'logiT | 

(ergs-1) - 1U \ 3 J 



(erg s 1 ) 



10 



53.1 / 

0.01 



(— 0.95±0.3O) 

-0.62±0.0.04) 
'(-1. 1410.10) 



BAT 



BATSE 



BATSE 



a Gold and silver sample from lukwatta et al. I <■-'<! Ill ' fitted ac- 
counting for the sample variance, 90% c.l. 
b 



Hakkil a et al.l on s ingle pulses 

Norris, Marani fc Bonnell j20Q0l) 



r 0.3-10 kcV 



(ergs- 1 ) 



_ ^q(50.82±0.20) x 



(-0.95±0.23) 



(9) 



where ii ag ,x is the rest frame peak lag calculated as the differ- 
ence in time between the flare peak times in the 0.3-1 keV 
and the 3-10 keV energy bands. The subscript x reminds 
that the time lag is calculated in the X-ray regime. We prop- 
erly accoun t for the sample v ariance following the method 
outlined by iD'Agostinil {2005). Errors are provided at 90% 
c.l. It is not possible to constrain the spectral peak energy 
and the spectral slopes of the majority of the flares of C10: 
for this reason we fitted the spectrum extracted around the 
peak time with an absorbed SPL within XSPEC and con- 
servatively quote the isotropic peak luminosity as obtained 
in the observed 0.3-10 keV bandpass (black bullets in Fig. 
[15)1 . In two cases, for GRB 060418 and GRB060904B the 
Band parameters were determined thanks to the BAT+XRT 
coverage (Table [2} : this allows us to plot the 1-10000 keV 
(rest frame) isotropic peak luminosity. However, the best fit 
relation was derived using their 0.3-10 keV values for ho- 
mogeneity (the two red bullets in Fig. [15] consistently lie 
above the expectation). In the cases of the GRB 060418 and 
GRB 090904B flares, the isotropic peak luminosity in the 
0.3-10 keV observed band underestimates the 1-10000 keV 
rest frame value of a factor 2-3. 

The lag- luminosity is one of the key relations which 
connects the GR B prompt temporal and spectra l proper- 
ties: discovered bv lNorris^ Marani fc Bonnell (2000) as a time 
integrated property of each particular burst, the relation 
was conclusively de monstrated to reflect pulses rather than 
bursts properties bv lHakkila et al] (2008). Figure [T5l shows a 
direct comparison between the flare and the prompt proper- 
ties in the lag-luminosity diagram: this is of particular inter- 
est since long and short burst s are known to occup y different 
regions of the plane (see e.g. iGehrels et al1l2006t l. Flares in 
long GRB are consistent with the long GRBs lag-luminosity 
relation. However, we should consider that: first, in Fig. 1151 
the lag of the prompt data is calculated using the Cross 
Correlation Function (CCF) to the entire BAT light-curve 
and consequently re flects a time integrated property (see 
lUkwatta et al.l l2010 for any detail). Second: for the prompt 
data the lag is defined as the time difference between light- 
curve structures in the 50-100 keV and 100-200 keV chan- 
nels. For the flares we calculated the peak lag between X-ray 
energy bands. In both cases the lag has been computed be- 
tween band-passes around the event spectral peak energy. 
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Rest frame Lag (s) 

Figure 15. Flare peak-lag peak-luminosity relation. Black filled circles: flares from CIO. The isotropic peak luminosity has been computed 
in the observed 0.3-10 keV energy band. The two flares (in GRB 060418 and GRB060904B) for which it was possible to estimate the 
Band parameters are marked with red bullets: in these cases the 1-10000 keV peak luminosity is plotted (see Table[2]|. In both cases the 
flare peak lag is reported. Red solid line: flare best fit peak-lag peak-luminosity relation: errors are provided at 90% c.l. Red dotted lines: 
best fit sample variance. Grey stars: prompt gamma-ray data corresponding to the gold and silver samples of long GRBs from Ukwatta 
et al., 2010. The isotropic peak luminosity has been calculated in the 1-10000 keV frame. The lag corresponds to the time difference 
between light-curve structures in the 50-100 keV and 100-200 keV channels. Grey dashed line: prompt best fit peak-lag peak-luminosity 
relation from Ukwatta et al., 2010. The shaded area marks the 68% region around this fit. Red dot-dashed line: best fit relation of the 
gold sample performed accounting for the sample variance Blue dashed line: Norris 2000 lag luminosity relation. 



Third: the prompt peak luminosity is calculated in the 1- 
10000 keV rest frame energy band, while for the flares the 
peak luminosity is calculated from the 0.3-10 keV (observed) 
bandpass (which is expected to be a factor of 2-3 below the 
1-10000 keV value). 

With these caveats in mind, it is remarkable that the 
bes t fit slope of the flare lag-luminosity is consistent with 
the lUkwatta et all (|2010h results based on BAT data (see 
Table [5]): L Pt i so oc t[~®^ 5±0 ' M ^ (red dot dashed line in Fig. 
I15p . It is not surprising that this relation overestimates the 
flare luminosity which is computed in a narrower energy 
window by a fa ctor ~ 5. Only a ma rginal consistency can 
be quoted with lHakkila et al.l (|2008h who reported an in- 
dex ~ 0.6. Our findings are inst ead fully consistent with the 
iNorris. Marani fc Bonnel |2000r i law (both in normalisation 
and index), with a power-law index 1.14 (blue dashed line in 
Fig. ll5[) : these authors reported lags between BATSE energy 
bands 100-300 kev and 25-50 keV. The same is true if we con - 
sider the lag- luminosity power-law index bv lSchaeferl (|2007T ) 
who reported a value of 1.01. The prompt lag-luminosity 
relations are summarised in Table [5] 



This result strongly suggests a common physical mech- 
anism producing both the GRB prompt emission and the 
X-ray flare emission hundreds of seconds later. The lag- 
luminosity relation has been proven to be a fundamental 
law extending ~ 5 decades in energy and ~ 5 decades in 
time. 

5.3 Flare parameters correlations 

At the 0th level of approximation, the flare phenomenology 
summarised in Table [3] reduces to a set of 4 independent re- 
lations (namely Ea. ll0irmil2l and ll3|) . This consistently lim- 
its the available parameters space for any theoretical model 
aiming at explaining the flare phenomena. 

The first fundamental relation expresses the concept 
that flares profiles tend to be self-similar both in energy 
and in time (C10): 

i r ~0.5i d (10) 

which implies: n ~ 4t2: the N05 profile contains a degree of 
freedom which is not required by the data, since the rising 
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and decaying time scales are not independent from one an- 
other; t r ~ w/3 ~ T2, td ~ 2u>/3 ~ 2r2 and ~ const. The 
pulse peak lag between two energy bands centered on _Bi and 
£ 2 reads: t Ug (E u E 2 ) = t s (£ 2 ) + 2t 2 ,e 2 - t s (£i) - 2r 2 , Bl ~ 
2(t2,b 2 — t 2 ,Ei) since t s {Eo) ~ t s (Ei) at the first level of 
approximation. T2 directly inherits the energy dependence 
of t r and t d (Table gj: t 2 ,e 2 = (E 2 /Ei)~ aT x T2,Bi where 
a T ~ 0.2 — 0.7. This implies ti ag ~ f(E\,E 2 ) x T2 where / is 
a factor of proportionality which depends on E\ and E 2 . A 
lag-width correlation is consequently expected since w oc r 2 
(or alternatively a ti ag — t r or ti ag — relations). 

The width-peak time linear correlation discussed by 
CIO, introduces the concept of temporal evolution, which 
is the key to interpret the GRB X-ray flare phenomenology: 

UI (X t pcak (11) 

which automatically implies t r oc t pca k and td oc t pca k; con- 
sidering the lag- width correlation of the previous paragraph, 
Eq. [TT] translates into a lag-tpeak correlation. 

The third fundamental relation links the flare peak-lags 
to their peak-luminosities, in strict analogy to the prompt 
emission: 

r ~, + -0.95±0.23 ft r,\ 

L p , iao oc t la& (12) 

Consequently, t r , td and w will be linked to the flare peak 
luminosity. In turn, the peak luminosity is expected to evolve 
with time to lower values, as found: L p ,i BO oc t~ e 2 a k although 
with a large scatter. 

The fourth relation again describes the temporal evolu- 
tion of the flare spectral properties: 

HR(t) oc t- Q ak (13) 

Flares become softer and softer as time proceeds. As a result, 
an inverse HR — t r , HR — td, HR — w and HR — ti ag relation is 
automatically built. This completes the set of 28 correlations 
(or lack thereof) of Tabled An important implication is that 
at the 0th level of approximation, the physics underlying the 
flare emission determines that out of 3 time scales describing 
the flare temporal profile (namely t t , td and t pea k ), only one 
degree of freedom survives. 

As a second level approximation, one should consider 
that t s = t s (E); f = f(Ei,E 2 ,a T ) which likely reflects 
the relative position of the spectral peak energy with re- 
spect to Ei and E 2 ; k — k(E) since the rise and decay 
times evolution with energy band slightly differs from one 
another(Fig. [6]). Any theoretical model aiming at explaining 
the flare emission is asked to be consistent with these set 
of findings, together with the flare spectral properties dis- 
cussed in the following paragraph. A critical revision of the 
existing models is in preparation. 

5.4 Temporal evolution of a flare spectrum 

5.4-1 E p (t) 

For the two flares with the highest statistics, we find that 
the flare spectrum is well described by an evolving Band 
function. In both cases the spectral peak energy evolves 
with time to lower values following an exponential decay 
which tracks the decay of the flare flux. After t pca k the N05 
profile is progressively dominated by the exp(— t/r 2 ) fac- 
tor, with the relative strength of the exp(— ri/t) dropping 



quickly. Remarkably, in both flares, the temporal evolution 
of the spectral peak energy and the flux seem to share the 
same e-folding time: for the flare detected in GRB 060418 
we have r 2 = 22.2 ± 0.4 s with r Ep = 20.0 ± 2.0 s; for the 
GRB 060904B flare we find r 2 = 29.5 ± 0.6 s in good agree- 
ment with te p = 29.0 ± 4.7 s. 

The evolution of E p (t) to lower values during the pulse 
decay time is one of the signatures of the prom pt emission 
as de monstrated by a number of studies (see e.g. iPeng et al.l 
120091 and references therein for a recent study on single 
prompt pulses). The spectral peak energy of prompt pulses is 
found to evolve following a power- law0S p (t) ~ t- s with an 
evolutionary slope S ~ 1 (see lPeng et al.|[2009l and references 
therein). The observed evolution of the flare E p (t) is instead 
much faster and inconsistent with the t _1 behaviour after 
rescaling the time to the beginning of the flare emission. Tak- 
ing the GRB 060904B flare which starts at t s ~ 124 s as an 
example, £ p (171s) ~ 4.7 keV would imply £ p (300s) ~ 1.2 
keV if the t _1 behaviour is assumed: the observations in- 
stead constrain the £ p (300s) to be well below the XRT 
bandpass, so that -E p (300s) < 0.3 keV. While this could 
be the result of a wrong choice of the power-law zero time 
(th e relevant time could be the ejection time instead of t s , 
see I Willingale et al.ll2010l for details), further investigations 
are required to establish if the faster E p (t) evolution is a 
peculiar feature of the flare emission when compared to the 
prompt phenomenology. 

5.4-2 Flares in the E p — E\ 

so and Ep — L p ^\ so planes 

The spectral analysis of Subsec. 14.31 allows us to constrain 
the properties of two X-ray flares in the E p — Ei so and 
E p — Lp.iso planes for the first time. Figure [T4l shows the re- 
sult: when compared to the spectral properties of the prompt 
emission of 83 G RBs with measured redshift and well con - 
strained spectra (|Nava et alj|2008l : [Ghirlanda et all 120091 ). 
the two flares show higher Ei BO and L Pt i so than expected 
starting from their rest frame spectral peak energy E p .i. If 
the flares are part of the prompt emission, then they are not 
expected to share the t ime integrated E v ,j — _Ei so best fit 
normalisation: however, iKrimm et al l l|2009l ) demonstrated 
that individual sequences of the same burst do follow the 
E p — Ei so relation with a higher normalisation (reflecting 
the fact that the total energy bud get is distributed over the 
sequence, see lKrimm et al]|2009l . their Fig. 14). The oppo- 
site is observed for the two flares. The E p — L p ,i so plane 
does not suffer from this effect, and time integrated prop- 
erties can be here directly compared to their time resolved 
counterparts. In this plane the two flares are barely inside 
the 3<r best fit area: however, the peak luminosity is poorly 
sampled below 5 x 10 50 erg s _1 and a deviation to lower spec- 
tral peak energies of the main sample cannot be excluded 
since the best fit relation is mainly established by high lu- 
minosity data points. In this case the X-ray flares would 
belong to the same E p — L Pl i SO as the prompt emission. Al- 
ternatively, and equally interesting, the flares do not follow 
the E Pl i — Lp.iso relation of the prompt data. This possibility 
is suggested by the fact that the two flares are outliers of 

6 The zero-time is chosen to be the starting point of the rising 
segment. 
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the E p 

solved spectral analysis is taken into account. Whether this 
depends on the chosen integration energy band (1-10 4 keV) 
is currently under investigation. 

For GRB 060418 the spectral pro perties during the 
prompt emission were measured (see iNava et al.l 120081 ; 
lAmati et al.ll200Sh : £jT ray = 572 ± 143 keV; E?-™ y = 
(13±3) x 10 52 erg; I/^ s r ay = 1.9 x 10 52 erg. The flare isotropic 
input comprises ~ 5% of the isotropic gamma ray prompt 
energy while the flare isotropic peak luminosity is ~ 3% the 
L^~^ y . The contribution of the flare emission is therefore 
well within the uncertainties estimated for the prompt pa- 
rameters: this results helps to understand how much the un- 
accounted flare emission contributes to the scatter affecting 
the E p a - E iBO relation (|Amati et al.ll200Sl ). 
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5.4-3 The flare E p (t) — L p ,is (t) relation 

The time resolved analysis of the two flares with redshift and 
BAT detection reveals that the rest frame spectral peak en- 
ergy E Pt i correlates with the isotropic luminosity L[ so within 
single flares. The same has been demonstrated to be true for 
the prompt pulses: the time integrated correlations are the 
result of the existence of similar time resolved correlations 
of the same parameters. Individual GRB pr ompt pulses are 
consis tent with the E p ,j — Lj ao as proved bv lGhirlanda et all 
l|2009h . and lOhno et all (|2009l ). Our work proves that 2 X- 
ray flares share this property and that flares in general are 
likely to be consistent with this behaviour. The best fit 
slopes of the two relations of Fig. [T3] are however steeper 
than the ~ 0.5 s l ope fo und for the prompt emission (see e.g. 
iGhirlanda et al J l2009h : whether this is something peculiar 
of the flare emission or not, needs to be understood with a 
larger sample. At the moment we note that the best fit time 
resolved E p ,i — Li so slope within single prompt pulses have 
been found to be different from the 0.5 val ue in some cases: 
see e. g. GRB 090323 and GRB 090328 in IGhirlanda et all 
(2009) their Fig. 3, where the E Pi i — L mo track seems to be 
steeper than 0.5. 

Intriguingly, the spectrum extracted at the onset of the 
flare in GRB 060904B is a clear outlier, being characterised 
by an E p ,j higher tha n expected: a similar result was found 
by lOhno et al.l (|2009h for the prompt emission. These au- 
thors performed a time resolved spectral analysis of the 
GRB 061007 prompt emission and concluded that the ini- 
tial rising phase of each pulse is an outlier of the E P: i — L[ so 
relation, with E P] i around twice the value expected from 
the spectral correlation. However, no difference between the 
rise and d ecay portions of t he pro mpt pulses has been re- 
ported by IGhirlanda et al.l |2009) from the analysis of 2 
Fermi GRBs. 

The presence of these spectral correlations during the 
flare emission tightly links the X-ray flare emission to the 
prompt phase; as a by-product, these results strengthen the 
interpretation of the spectral energy co rrelations as man ifes- 
tation of the physics of the GRBs (see Na va et al. I l2008l and 
references therein for a detailed discussion of this topic). 



Figure 16. Red line: classification of the GRBs belonging to the 
flare sample of C10 according to their X-ray afterglow morphology 
(see Evans et al., 2009 for details). Black line: GRBs belonging to 
the Evans et al. 2009 sample detected in the same period of time. 
This distribution has been re-normalised for the total number of 
GRBs of the flare sample. 



5.5 Flares and light-curve morphology 

The entire set of findings presented in this work establishes 
a strong link between early X-ray flares and prompt pulses. 
Flares are consistent to be independent episodes of emission 
superimposed to the contemporary steep decay or afterglow 
components: while the first has been proven to be gener- 
ated by the t ail of the previous puls es dKumar fc Panaitescul 
2000; see e.g. IWillingale et al. 2010) and is therefore related 
to the physical mechanism which gives rise to the prompt 
emission, the lat ter is likely to be a co mpletely independent 
component (e.g. iMargutti et al.ll2010l ). 

X-ray afterglows with no breaks are likely to be dom- 
inated by the afterglow c omponent from th e beginning of 
the XRT observations (see lLiang et al.| [2009 for a dedicated 
study). We therefore expect the SPL X-ray afterglows to 
show a limited flaring component if flares really constitute 
an independent component (flares would be perhaps present 
but hidden by the contemporaneous afterglow emission). 
Figure [16] confirms this expectation: the no-break class is 
clearly under-represented in the flare sample, while the vast 
majority of the detected early time flares reside in one-break 
or canonical X-ray afterglows. The probability that the two 
samples are drawn from the same population is evaluated 
to be as low as 0.9% using a KS test. If this is due to 
the brightness of the afterglow emission, than the conclu- 
sion is that flares have an independent origin: in the stan- 
dard scenario this means that flares are not produced by 
the external shock. Another equally interesting possibility 
is that flaring emission is quenched (instead of being hid- 
den) in SPL-afterglow GRBs. Distinguishing between these 
two possibilities is beyond the scope of this work and will 
be the subject of a forthcoming investigation. 
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Figure 17. Flare paradigm. Dashed curves: high energy profile. 
Solid line: low energy profile. High energy profiles rise faster, de- 
cay faster and peak before the low energy emission. As time pro- 
ceeds, flares becomes wider, with higher peak lag, lower peak 
luminosities and softer emission. However, at the Oth level of ap- 
proximation the self-similarity is preserved both in energy and 
in time (with t r /t^ ~ 0.5). We stress that this is an average 
behaviour. The thick solid arrow underlines the dependence of 
the flare properties on the time, in strong opposition to prompt 
pulses. 



the flare the lower is the rate of evolution from one energy 
band to the other; 

• The flare spectral peak energy E p (t) evolves to lower 
values following an exponential decay which tracks the de- 
cay of the flare flux. The detected evolution is faster and 
inconsistent with the ~ t _1 behaviour even when the zero 
time is re-set to the beginning of the flare emission; 

• The two flares with best statistics show higher than 
expected Ei so and Lp,i so values when compared to the E Pi i — 
.Eiso and E Pi i — Li so prompt correlations; 

• The rest frame spectral peak energy _E P; i correlates with 
the isotropic luminosity Li so within single flares, giving rise 
to a time resolved E Pl i — L[ so correlation; 

• The flare emission preferentially builds up at lower en- 
ergies: flares do not seem to be consistent with the Pulse 
Start Conjecture; 

• Among the different types of X-ray afterglow light- 
curves, the simple power-law afterglows are under- 
represented in the flare sample. Flares are preferentially de- 
tected superimposed to one-break or canonical light-curves. 

The strict analogy between the prompt pulses and flares 
phenomenology strongly suggests a common origin of the 
two phenomena. 



6 SUMMARY AND CONCLUSIONS 

The GRB X-ray flare phenomenology shows a list of prop- 
ertie s strictly analogous to the gamma-ray prompt emission 
(see lHakkila et alj [2008 for a recent study on individual 
prompt pulses). However, it differs in one crucial aspect: 
flares evolve with timfl Flares evolve with time to lower 
peak intensities, larger widths, larger lags and softer emis- 
sion. Fig. [17] illustrates the flare paradigm. As the prompt 
emission pulses, flares have correlated properties: short-lag 
flares have shorter duration, are more luminous and harder 
than long-lag flares. In particular the following properties 
add to the lists of Sec. \T\ 

• Flares define a lag-luminosity relation: 

L°t" 10kcV oc ir°- 95±0 - 23 . The best fit slope is remarkably 
consistent with the prompt findings (see lUkwatta et al.ll2010l 
and references therein for a recent study); 

• The lag is found to correlate with the flares width while 
it is inversely correlated to the flare spectral hardness; 

• The flares temporal profiles in different band-passes are 
only a nearly exact time stretched version of one another: the 
rise and decay times evolve following slightly different power- 
laws in energy: t r oc E~ atr , td oc E~ atd with cttr > cttd- The 
result is that flares are on average more asymmetric at high 
energy; 

• The rate of evolution of a flare profile in different band- 
passes anti-correlates with its spectral hardness: the harder 

7 We parenthetically note that this difference could be due to 
the fact that the relevant time for a peak (prompt pulse or flare) 
is the time since the shell ejection in a standard internal shock 
scenario, and not the time since trigger. While these two time 
scales are completely unrelated for the prompt pulses, they are 
more closely related for flares (since the ti me is large the dif f erence 
between the two becomes negligible). See [Willingalc et al.l i2010T ) 
and Willingale et al., in prep, for further details. 
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Table Al. Flares best fit parameters according to Eq.[T]and derived quantities (Eq.|2l[3]|4]and[5jl. The reported uncertainties are calculated 
accounting for the entire covariance matrix for each fit. From left to right: GRB name, energy band of investigation, redshift, normalisation, 
start time, first shape parameter, second shape parameter, peak time, width, asymmetry, rise time, decay time, chi squared, degrees of 
freedom. In each column, a -1 indicates the absence of the value in the literature (for the redshift) or the impossibility to constrain the 
value from the fit. A negative error indicates that the parameter was first left free to vary in the fit and then frozen to the reported value 
to constrain the errors associated to the other fit parameters. 
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0.291 
0.290 



a 0.3-10 keV 
b 0.3-1 keV 
c 1-2 keV 
d 2-3 keV 
e 3-10 keV 

/ First flare detected by the XRT. 
9 Second flare detected by the XRT. 
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Table A2. Best fit parameters derived from the spectral modelling of the XRT and BAT data of GRB060904B. A photo-electrically 
absorbed models (tbabs*ztbabs within xspec). Two different spectral models have been used: a simple power-law (SPL) and a band 
function with the peak of the vF v spectrum as free parameter. From left to right: name of the interval of the extraction of the spectrum: 
XRT+BAT stands for a joint BAT-XRT data fitting; start and stop time of extraction of each spectrum; model used; intrinsic neutral 
Hydrogen column density; best fit low and high energy photon indices for a Band function or best fit photon index Y for a SPL model; 
spectral peak energy; isotropic luminosity; statistical information about the fit. The * symbol indicates an apparent trend in the residuals 
of the fit. A joint fit of all the spectra giv es JVh,z = 0.5 ± 0.1, x 2 /dof=579. 76/567, P-value=35%.The Galactic absorption has been frozen 
to 1.13 X 10 21 cm -2 l|KaIberIa et al.ll2005l l. Errors are provided at 90% c.l. 
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< 1.01 


0.17 + 0.02 


89.96/78 
85.34/78 


0.17 
0.27 


7 


XRT 


224 


249 


SPL 
Band 


0.88 + 0.12 
0.52 


1.00 


3.05 + 0.13 
3.05 ± 0.22 


< 0.4 


0.087 + 0.003 


79.94/78 
91.89/77 


0.42 
0.12 


8 


XRT 


249 


320 


SPL 
SPL 


0.71 + 0.09 
0.52 




3.60 + 0.18 
3.28 ± 0.08 




0.030 ± 0.003 


92.80/70 
105.94/71 


0.04 
0.01^ 



a The presen c e of a spectral feature in excess of the SPL component in this interval of time is discussed in iMargutti et al . (2008) and 
iMoretti et all J2008h . 

Table A3. Same as Table \M\ for GRB 060418. A joint fit of all the spectra give s AT H , Z = 0.6 + 0.1, x 2 /dof=318.43/327, P-value=62%. The 
Galactic absorption has been frozen to 8.81 x 10 20 cm -2 l|Kalberla et alj|2005t l. Errors are provided at 90% c.l. 





Interval 


U 
(b) 


*/ 
(b) 


Model 


(10 22 cm" 2 ) 


"B 


Mr) 


-^peak 

(keV) 


Lisa 

(10 50 erg s" 1 ) 


X 2 /dof 


P-value 


1 


XRT+BAT 


122 


131 


SPL 
Band 


2.2 + 0.50 
0.59 


0.90 


2.0 + 0.1 
2.37 ± 0.22 


c i+Q.9 
°-l-0.7 


4.2 + 0.4 


126.03/64 
74.26/64 


6 x 10~ 6 
0.18 


2 


XRT+BAT 


131 


140 


SPL 
Band 


2.6 + 0.4 
0.59 


0.90 


2.1 + 0.1 
2.33 + 0.13 


4 =7+0.5 
4 -°'-0.6 


6.6 + 0.6 


143.20/103 
86.16/103 


5 x 10~ 3 
0.88 


3 


XRT+BAT 


140 


148 


SPL 
Band 


1.91 + 0.42 
0.59 


1.00 


2.29 ± 0.12 
2.50 + 0.40 


9 9 +0.3 


2.6 + 0.2 


62.85/58 
62.30/58 


0.31 
0.33 


4 


XRT+BAT 


148 


170 


SPL 
Band 


1.18 + 0.17 
0.59 


1.00 


2.44 + 0.09 
2.44 ± 0.11 




1.0 + 0.1 


96.44/107 
96.32/107 


0.76 
0.76 



E Pli (keV) 




